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Abstract 

We investigate the electronic transport through a single molecule in a strong electron-phonon cou- 
pling regime. Based on a particle-hole transformation which is made suitable for non-equilibrium 
situation, we treat the pair tunneling and cotunneling on an equal footing. We propose an experi- 
mental setup to enhance the visibility of pair tunneling, which has no Franck-Condon suppression. 
We also discuss the shot noise characteristics. 
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I. INTRODUCTION 



Many transport features of single molecules are ascribed to the coupling between elec- 
tronic and vibrational modes 1,2 . Especially, the electron-phonon coupling makes the effective 
charging energy U reduced and even negative. The negative U makes the double occupancy 
of the molecular level energetically favorable, and allows the tunneling of pairs. Although 
most studies on the negative U have focused on the Kondo regime^ 1 ^, a recent work- pointed 
out that the pair tunneling can significantly affect transport above the Kondo temperature 
(Tk) as well, and leads to several very unusual features^. In particular, the resonance width 
of the pair tunneling conductance was only determined either by the temperature (linear) 
or by the bias voltage (non-linear); but not by the lead-molecule hybridization as in normal 
sequential tunneling of single particles. Further, the pair tunneling conductance was exactly 
twice the normal cotunneling contribution. 

However, we note that there are several physically important questions that have not been 
addressed yet. Firstly, there is no justification why the pair tunneling be treated separately 
in a sequential tunneling picture apart from the normal cotunneling processes. In this work, 
we use a different approach and treat both on an equal footing. Secondly, it is a general 
observation that the pair tunneling is usually subject to the exponential Franck-Condon sup- 
pression. Then, how can we observe the pair tunneling processes experimentally in a broad 
background of normal cotunneling contributions? Here we propose an experimental setup 
enabling one to observe the pair tunneling. Thirdly, is there a simple physical explanation 
why the pair tunneling conductance is exactly twice the normal cotunneling contribution? 
Can one attribute it to the double charge of the pair (like in a tightly bound pair objects)? 
We find that this is not the case. 

In addition to the previous questions on average current (conductance), for a better 
understanding of the nature of the pair tunneling, one can also investigate the fluctuations 
of the current. In many systems, the effective charge of the elementary excitations has been 
identified by the shot noise characteristics^^. Interesting questions would then be: Should 
the pair tunneling events give the fluctuations corresponding to the Fano factor equal to 2? 
How different are the pairs tunneling through the device from tightly bound pair objects? 
In this work, we address all these questions and provide further physical insights to the 
pair-tunneling transport (above Tk), based on particle- hole transformation. 
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The paper is organized as following: In Section [TT] we establish the particle-hole trans- 
formation, which enables us to treat the pair-tunneling processes on the same footing as 
the normal single-electron cotunneling processes. This transformation, augmented with 
the exchange of the electrode indices, is suitable for non-equilibrium (as well as equilib- 
rium) transport. In Section II I II we evaluate all the relevant cotunneling rates based on 
the particle-hole transformation. This will make clear the difference between the underly- 
ing physical processes governing the pair tunneling and normal single electron cotunneling. 
Moreover, we clearly point out how the pair tunneling transport is exponentially suppressed 
by the Franck-Condon effect in realistic molecular devices, which was overlooked in the pre- 
vious work. 6 Section IIVI is then devoted to reproduce the results on average current (and 
conductance) of previous work to demonstrate the efficiency of our method. The current 
fluctuation noise and investigate the contribution to it from the pair tunneling processes in 
Section |VJ In Section IVII a realistic experiment is proposed where one can observe the pair 
tunneling physics without suffering from the Franck-Condon suppression, and provide the 
detailed analysis of the setup. Section IVHl concludes the paper. 

II. PARTICLE-HOLE TRANSFORMATION 

We describe the molecular device with an Anderson-Holstein model (AHM), whose Hamil- 
tonian has the form H = + Hr + Hm + Hj>- The left (L) and right (R) electrodes are 
described by the non-interacting electrons 

H e = Y,^k-eVe)c\ klT c tka (£ = L,R) (1) 

kcr 

Here c\ ka creates an electron with momentum k and spin a on the lead I = {L,R}. The 
molecule has both the electronic degrees of freedom, described by the fermion operators d 
and d\ (n a = d\d a , h = n-j- + hi), and the vibrational mode with frequency uo, described by 
the boson operators a and a): 

Hm = cadi-dp + U h^hi — Xhuo(h — l)(d^ + d) + TiuqcJcl , (2) 

where A is the dimensionless electron-phonon coupling constant between the two degrees of 
freedom. The electron tunneling between the each lead and the molecule is expressed as 

Ht = E ( T ft4A + h.c) . (3) 

Ika 
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The molecule-electrode couplings are characterized by the hybridization parameters Tg = 
2n J2k \Tek\ 2 5(E — e ka ). As usual, we ignore the weak energy dependence of IY We also put 

r = r L + r R . 

The molecular part H M is diagonalized as 



H 



M 



,+S 



e s + hujncJa 



(4) 



by the canonical transformation S = X(n — l)(a* — a).-^ The renormalized molecular 
level €d and on-site interaction U are given by = eo + X 2 huo and U = Uq — 2X 2 hujQ, 
respectively. In the strong coupling limit (A 2 ^> Uq/2Tiujq) the on-site interaction U becomes 
negative so that a double occupancy of the molecular level is energetically favored to a 
single occupancy. Indeed, by removing the phonon degrees of freedom with a Schrieffer- 
Wolf (SW) transformation, it was shown that electrons are allowed to tunnel in pairs into 
the molecule^; see Fig. [1] (a). However, it should be stressed that the situation considered 
in Ref. f] is rather special and in fact equivalent to a negative-f/ Anderson model {isotropic 



Kondo model). Earlier, it was shown in Ref. 
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that the model is, in general, equivalent to 



an anisotropic Kondo model (not equivalent to a negative-t/ Anderson model).™ We will 
come back to this issue later. 

It has been proved useful to map a negative-f/ impurity model to an equivalent model 
with positive interaction by a particle-hole transformation (PHT)^ 2 -. In the same spirit, we 
directly apply the PHT to the AHM keeping the phonon modes^. We first choose an one- 
to-one correspondence k i— > k such that t~ k = — (we assume symmetric conduction bands). 
Following Refs. , we then make the PHT for down spins 



d i ^ d\ , c Lkl 



2kl ' 



c Rkl 



5 f - 
lki ' 



(5) 



keeping the up spins unchanged 



d T h-> d T , c Lkr 



(6) 



It is emphasized that the lead indices for the down spins have been exchanged (L — > 2, R — > 
1) compared with those for the up spins [L — > 1,R — > 2). This is a small yet important 



difference between the mapping here and that in Refs. 
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12c the lead-index exchange is not 



necessary at equilibrium. Under this transformation, Hl + Hr is transformed to Hi + H2 
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with 

Hi = ^(efca - eVe a )cl ka c ikf7 , (7) 

ka 

and H T to 

H T = J2 {^ikAkA + h.c) . (8) 

Ika 

The effective bias voltages eV^ become, in general, spin-dependent: eVif = eV^, eVij_ = 
—eVu, eV 2 i = eVu eV^j, = — eVj,. For the symmetric bias, which is assumed here, they are 
spin-independent: eVy\ = eVii = eVb/2, eV 2 ^ = eV 2 ± = —eV^/2. The tunneling amplitudes 
Ttko are also spin-dependent and given by 

T. rri rri rri* rri rji rri rri* / f\\ 

l^ — J-LkU J-lkl — -l-Flkl ' J2fc T — J RfcT > l2k l = ~ 1 Lkl- ^ ' 

Their spin-dependence disappears only for symmetric junctions (Tl = Th = F/2): = T/2 
[I = 1, 2). Later, we will see that the junction asymmetry affects significantly the transport 
properties of the device. Finally, the molecular part now takes the form 

(10) 

with S = A(n^ — ni)(a'— a). The effective on-site interaction U = —U thus becomes positive, 
and the molecular level Ida = U/2 + cxAz/2 has an effective Zeeman splitting 

A z = 2e d + U = 2e + U . (11) 

The eigenstates are given by |0)|m) (with energy E 0tm = mhuo), \2)\m) {E 2j7n = mhujo), 
and D(a\)\a)\m) {E a>m = mfiuio + Ida), where |0) (empty), |2) (doubly occupied), and \a) 
(singly occupied by spin a) are electronic states of the molecule, \m) is the phonon state, 
and D(a) = exp [a(a! - a)] 

In the strong coupling regime (A 2 ^> U /2huj ), the equivalent model has U ^> T and 
Qct — T < 0; the so-called local-moment regime. The transport thus occurs only through 
cotunneling processes (above TV). The PHT maps the electronic state |0) of the molecule in 
the original model to | J) = d\\0) in the equivalent model; likewise, ||) i— > |2), | j) i— > |0), and 
1 2) l— > IT)- The usual cotunneling process in the original model thus corresponds to the spin- 
preserving cotunneling (SPC) in the equivalent model, and the pair tunneling [Fig. 12(a)] 
to the spin-flip cotunneling (SFC) [Fig. [2(b)]. Namely, unlike in Ref. Q, in our picture 



edad\d a + Un^ni + a) a 
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FIG. 1: (color on-line) An illustration of the particle-hole transformation. A pair tunneling in the 
original model (a) corresponds to a spin- flip cotunneling in the equivalent model (b). To emphasize 
the exchange of the electrodes in the particle-hole transformation, different widths (colors) have 
been used for electrode bands. Spin-preserving cotunneling processes (not shown) remain the same 
in both models. In (a), shown on the molecule site is the single-particle energy level for the double 
occupancy. 

all the relevant processes are treated on an equal footing, in terms only of "cotunneling". 
This will allow us to infer further insight into the pair-tunneling transport based on the 
well-established theory of cotunneling transport 15 . 

A few further remarks are in order: (i) In most experiments, the two leads are identical. 
But the mapping applies to general cases, including ferromagnetic leads^S. (ii) The trans- 
formed Hamiltonian arises from a mathematical mapping. Yet, the model itself is physical 
(experimentally realizable). Among others, an interesting realization will be the capacitively 
coupled double quantum dot (CCDQD); see Fig. HI 
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FIG. 2: Cotunneling rates in the equivalent model, (a) Spin-preserving cotunneling, corresponding 
to the usual single-particle cotunneling in the original model, (b) Spin-flip cotunneling, corre- 
sponding to the pair tunneling in the original model. 

III. COTUNNELING RATES 

Both the average current and the shot noise characteristics can be obtained from the rate 
equation 



+ J2w^P a/ (N-l)+J2 W ^'MN + l). (12) 



Here P (T (N,t) is the probability to find the molecule in the electronic state \a) and N 
electrons in the lead 2 at time t. W^, is the the rate of the cotunneling process where an 
electron is transferred from lead £' to I and, at the same time, the molecular state changes 
from | a') to \a). The cotunneling rates for the spin-preserving and spin- flip processes are 
given respectively by 



F{eV ia - eV t , a ) 



1 1 



X 



Fjida ~ (-da + eVto- - eV e 
2nTi 



(13) 



(14) 



where u = —U/2Hujq and F(E) = E/(e l3E — 1). One can immediately notice the differences 



between the two rates (see also Ref. l25l ) . (i) As illustrated diagrammaticaliy in Fig. [2],— 
the SPC rate is an incoherent (classical) sum of two contributions (because the spins of 
the conduction electrons are different between the two corresponding diagrams) while the 
SFC rate is a coherent sum of two indistinguishable contributions. The latter gives rise to 
a constructive interference for the SFC rate, and hence becomes exactly twice W%%. at 



eVb ~ and Az ~ 0. This explains why the pair-tunneling conductance is twice the normal 
cotunneling contribution, (ii) They have different Franck-Condon (FC) factors 

*•*(.., A) s< r*f; (±A 7" I (15) 

/ — ' ml m + u 

m=0 

The FC suppression arises at two different levels. According to the FC principle^, each 
virtual tunneling amplitude itself has the FC suppression factor: (a, 0\Ht\qz, m) ~ e~ A2 / 2 
and (a, m\H'r\o', 0) ~ e _A2//2 (a = 0, 2). For the overall FC factor, however, the contributions 
from various intermediate states should be all summed up. (i) When Huj ^> —U/2 (u <C 1), 
\a,m = 0) is the only intermediate state contributing to the overall cotunneling amplitude. 
In this case, the overall FC factor is e~ A2 for both W^ c and W^. This is the limit considered 
in Ref. [f], corresponding to the isotropic Kondo model (negative-U Anderson model), (ii) 
However, if hu>o <C —U /2 (u ^> 1), all the intermediate states \a, m) with higher vibrational 
energies give finite contributions. The sum of these contributions just amounts to cancel the 
FC suppression factor from the individual virtual tunneling for W^, while it is not the case 
for W^. Consequently, the SFC (the pair tunneling) is exponentially suppressed compared 
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and 



with the SPC. This is the case examined by most authors, including Refs. 

We define the relative FC factor by 7 = T-{u, A)/JF + (m, A). In typical experiments with 
Ceo molecule^, Uq ~ 300 meV and Kuoq ~ 5 meV. Therefore, the condition for the case 
(i) above is hardly satisfied, and 7 can be significantly (even though not exponentially) 
smaller than 1, and it may be difficult to observe the pair tunneling. Below (Section IVip we 
will propose an experiment where the pair tunneling is not subject to the Franck-Condon 
suppression. 



IV. CURRENT 

The average dc current is given by I = J2 a (T da ) with T da = (T R(7 -T La )/2 = (I 2a -ha)/2, 
where 1^ = edN^/dt and Ie a = edNg a jdt. Referring to the PHT (FigHJ), one can identify 
the two contributions to the total current I = e(J c + J p ), namely, 

J c = (J TT + J u )/2 (16) 

from the usual cotunneling and 

J p = JnP, + J iT P T (17) 
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from the pair tunneling. Here J a 



W. 



21 



W%, and P a 



J2 N P a {N,t = oo) is the 



stationary probability distribution of the molecular state. Putting the cotunneling rates 
in Eq. ffT3l) . one can easily reproduce all the results in Ref. [fi, in particular, the linear 
conductance (up to the FC factor 7 2 ) 



G 



2e 2 r L r R 
h 



j 2 U 2 



P(2e d + U) 



_ e 2( ed + [/)2 2sinh[/?(2e d + [/)] 



+ 



l-f(2e d + U) f(2e d + U) 

4 (^ + uy 



(18) 



Here we just clarify the questions raised at the beginning. In the original sequential-tunneling 
treatment of the pair tunneling in Ref. O, it is not clear why and to what extent the coherence 
between the subsequent pair-tunneling events can be ignored. In the present picture, the 
pair tunneling and normal cotunneling are treated on an equal footing, all in the cotunneling 
picture. Thus, in order to go beyond the sequential-tunneling treatment of the pair tunneling, 
viz., to the purely coherent resonant tunneling of pairs, which leads to the Kondo effect^, 
one has only to go to higher orders. Further, pair-tunneling contribution is exactly twice 
(for 7 = 1) the usual cotunneling contribution due to the interference in the SFC process 
(see above and Fig. [2]). 



V. NOISE 



The current noise spectral density is given by (at sufficiently low frequencies uo <C —U /2) 

/>oo 

S(u) = / dr e + ^(ST d (r)5I d (0) + 5T d (0)5I d (r)) (19) 



where SI d = I d — (I d ). The Fano factor, defined by S(0)/2eI, is a representative characteristic 
of the shot noise (eV& 2> fcsT). It may reveal not only the correlated transport but also the 
effective charge of the carriers^"^. Using the quantum regression theorem^yiL 2 ^ anc j t ne 
rate equation (fl2"]l . we obtain current-current correlation function (/c#T eVb —U/2 and 
A z < -U/2) 



S(lu) = 2el - 4e 2 
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1 + u 2 /W t 2 



p J 



(20) 



where W p = Y^Wa The Fano factor thus has the form 



2el 



Jn z J 



ir 



Jn ~\~ Jr) 



(21) 



with J c and J p given by Eqs. (TTB1) and ([T 




FIG. 3: (color on-line) (a) Fano factor as a function of A z /eVb for different values of a from 
(black solid line) to 0.9 (violet empty square) in steps of 0.1. a = (Tl — Tn)/T is the junction 
asymmetry parameter. Here eVj, = 500/csT and 7 = 1. (b) Fano factor with respect to I p . 

The shot noise characteristic is summarized in Fig. [31 It shows several interesting features: 
(i) The deviation from the Poissonian shot noise [S(0)/2eI = 1] comes entirely from the pair 
tunneling. This can be easily understood since the SPC in the equivalent model is an elastic 
process, which is known to be Poissonian^. (ii) For the symmetric junctions (Tl = Tr), 
the Fano factor is exactly 1 at resonance (Az = 0), and increases rapidly with \Az\ (up to 
\Az\ ~ eVb). The usual sequential tunneling at resonance across symmetric junctions gives 
a Fano factor of 0.5^. One may be tempted to interpret it as 0.5 x 2 with the factor of 2 
for the effective charge 2e of pairs. However, this cannot be justified since the Fano factor 
is relative to the total current / = e(J c + J p ), but not to the pair-tunneling contribution 
I p = eJ p . For definiteness, let us examine the pair Fano factor S(0)/2eI p , as shown in Fig. [3] 
(b). It shows that the pair tunneling in the system cannot be interpreted as a tunneling of 
tightly bound pair objects.— Instead, we interpret it again in the cotunneling picture of the 
equivalent model. The SFC is an inelastic process. At Az = 0, the two channels for the 
inelastic cotunneling are equal, and give no additional fluctuations (the Fano factor of 1). 
As Az increases, one channel carries larger current than the other, which gives additional 
fluctuations^, (iii) The shot noise turns out to be very sensitive to the junction asymmetry. 
Unlike single-particle sequential tunneling, the pair-tunneling even gives sub-Poissonian shot 
noise. This is also different from the inelastic cotunneling noise in usual QDs— , but this 
difference comes from the exchange of the lead indices in the PHT. 
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FIG. 4: (color on-line) A schematic of capacitively-coupled double quantum dots. The majority 
carriers through QD 1 (upper) are the electrons (red filled circle) and those through QD 2 (lower) 
are the holes (red empty circle). The gate voltages coupled to QDs are not shown for simplicity. 



VI. TUNNELING OF ELECTRON-HOLE PAIRS 



Now we propose a experimentally feasible set-up to observe the pair tunneling transport. 
It consists of capacitively coupled double quantum dots, where each dot is connected to its 
own leads via tunneling junctions and two dots are coupled in parallel capacitively; see Fig. HI 
The coupling capacitance between two dots is denoted by C\, to be distinguished from the 
self-capacitance Cq of each quantum dot(QD), which includes the junction capacitance and 
the gate capacitance. The tunnel junctions are specified by the hybridization parameter 
Fe={L,R},i={i,2} between QD i and lead I. Since capacitance is relatively insensitive to the 
sample fabrication geometry (unlike exponentially sensitive hybridization), for simplicity we 
assume that the capacitances are symmetric over the upper and lower branches and for the 
left and right junctions. The gate voltages are applied oppositely to the two QDs, so that 
the gate-induced charges are opposite to each other: 

n g ,i = -n a , 2 = n g . (22) 

The effective carriers, say, through QD 1 are then electrons while those through QD 2 are 
holes, provided that the QDs are made of semiconductors^. When the QDs are ultra-small 
metallic grains^2£, the "electrons" here are the excess electrons and the holes are the deficit 
electrons with respect to the mean background charge. 
The electrostatic energy of the double QD is given by 

E c (ni, n 2 ) = £i(m + n 2 - 2n g f + E Q (m - n 2 f (23) 
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FIG. 5: The Coulomb energy ratio E\/Eq versus the capacitance ratio C\/Cq. Note that for 
arbitrary C\/Cq, Eq is always bigger than Ei, leading to the negative interaction energy Uc- 

where n\ is the number of excess electrons on QD 1 while n 2 is the number of excess holes 
(the number of deficit electrons) on QD 2. In Eq. fl23|) . the two Coulomb energy scales E 
and Ei have been defined by 

E = (24) 

and 



E x = 



8Ci 1 + C /2d 
We note that the ratio Ei/E is given by 

Ei __ _ Co/2Ci 
£"0 



(C /2Ci) - (Co/2d 



(25) 



(26) 



1 + 

and that i?i is always smaller than Eq for arbitrary ratio C\ / Co as illustrated ed in Fig. [51 
In particular, for C\ ^> Co, E\/Eq becomes vanishingly small. 

In this large coupling regime, Eq. (|23|) shows that the electrostatic energy is very low 
when there are same number of electrons and holes in QD1 and QD2, respectively, that is, 
when ni =712- 111 other words, it is energetically favorable for this system to form electron- 
hole pairs in two dots. Indeed, in the Coulomb blockade regime (E , E x ^> T, k B T), if we set 
the gate voltage, so that n g = 1/2, then one can identify four lowest-energy charge states of 
the DQD, \ni,n 2 ) = |0,0), |1, 1), |1,0) and |0, 1), whose energies are given by 



E c (0, 0) = E C {1, l) = E 1 , E c (l, 0) = E c (0, 1) = E 



(27) 



Therefore, Uc = 2(E\ — Eq) < plays the role of the negative interaction energy U of the 
molecule and e c = (E — E x ) > is the counterpart of the energy e d of the singly-occupied 
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states in the molecule. We emphasize that for arbitrary Co/C\, one always has Uc < and 
ec > 0; see Eq. ff26l) and Fig. [51 Of course, in order for the electron-hole pair state |1, 1) 
(together with the empty state 1 , 0) ) to dominate the other charge states in the transport, 



it is required that 



which in turn gives the condition 



|£He c >r, (28) 



c] ^ 2(Eo/T - 1) • (29) 
In typical experiments on small metal-grain QDsM^^-, Eq/T > 5. It then suffices that 

I » To ' < 3 °> 

Large capacitive couplings have already been realized experimentally. C\/Cq ~ 5 in 
nanoscale metallic grains^, fully satisfying the condition (!30l) . On semiconducting quantum 
dots^i Ci/Cq ~ 0.28(0.34), satisfying the condition fl30l only marginally and hence reducing 
the pair-tunneling contributions slightly from the ideal values. A very recent measurement 
of the cross-correlation noise through CCDQD^ 2 -, even though the coupling capacitance in 
this experiment was rather small {Ci/Cq ~ 0.1) and pair-tunneling contribution cannot be 
dominant, suggests that further detailed studies of the CCDQD are worthwhile. 

Given the condition fl30|) satisfied, the calculations of the current and noise are exactly 
the same as in the molecular case with only one exception. The exceptional difference 
comes from the fact that, for example, Fl,i and Tl,2 can be significantly different in realistic 
experiments. These junction anisotropies change, e.g., the conductance in Eq. (1T81 to 



q _ 2e 2 T Ly iT R 2 + r L 2 r_R,i 



h 



U 2 (3{2e d + U) 

e 2 d {e d + Uy2smh[(3(2e d + U)} 



2e 2 T L iT R i + T L 2 T R 2 
h 2 



l-f(2e d + U) f(2e d + U) 
4 (e d + U) 2 



(31) 



Note that the FC factor 7 2 does not appear here, and hence for reasonably symmetric 
junctions with (r^r^ + T^r^i) w (Xl,i^r,i + r^^r^^) the pair-tunneling contribution 
can be clearly seen on top of the broad normal-cotunneling background. The junction 
anisotropy is expected to enhance the current fluctuations and hence the super-Poissonian 
nature of the noise characteristics. 
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VII. CONCLUSION 



By exploiting the particle-hole transformation, we studied the pair tunneling through a 
single-molecule transistor on an equal footing to the normal cotunneling. We have clarified 
the nature of the pair tunneling by revealing new features of the shot noise characteristics 
as well as reinvestigating the average current. We also respected the general observation 
that pair tunneling is subject to stronger FC suppression than the normal cotunneling, and 
proposed an experimental setup to enhance a visibility of the pair tunneling. 
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